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Abstract. We investigate nonlinear dynamics near an unstable 
constant equilibrium in the classical Keller-Segel model. Given any 
general perturbation of magnitude 5, we prove that its nonlinear 
evolution is dominated by the corresponding linear dynamics along 
a fixed finite number of fastest growing modes, over a time period 
of In | . Our result can be interpreted as a rigourous mathemati- 
cal characterization for early pattern formation in the Keller-Segel 
model. 



I. Growing Modes in the Keller-Segel Model 

The goal of this section is to review the well-known instability crite- 
rion for the classical Keller-Segel model, which describes directed move- 
ment of microorganisms and cells stimulated by the chemical which 
they produce themselves. The Keller-Segel system takes the form 

(1.1) [/ t = -V(-/iW + x^W), 
V t = V (DVV) +fU- kV, 

where U(x,t) is the cell density, V(x,t) the chemo-attractant, fi > 
the amoeboid motility, \ > the chemotactic sensitivity, D > 
the diffusion rate of cAMP, / > the rate of cAMP secretion per 
unit density of amoebae, k > the rate of degradation of cAMP in 
environment. 

We assume Neumann boundary conditions for U(x,t) and V(x,t), 
in a d- dimensional box x G T d = (0, 7r) d , d = 1, 2, 3, i.e., 

dU dV 

(1.2) — = — = 0, at Xi = 0, 7r, for 1 < % < d. 

OXi OXi 

A uniform constant solution 

U(x,t) = U, V(x,t) = V 
forms a homogeneous steady state provided 



(1.3) 



fU = kV. 
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In this article, we study the nonlinear evolution of a perturbation 
u(x,t) = U(x,t) - U, v(x,t) = V(x,t) - V 
around [U, V], which satisfies the equivalent Keller-Segel system: 

(1.4) u t = /iV 2 w - xUV 2 v - yV(wVv), 

(1.5) vt = DV 2 v + fu- kv. 

The corresponding linearized Keller-Segel system then takes the form 

(1.6) u t = fiV 2 u - xUV 2 v, 

(1.7) v t = DV 2 v + fu- kv. 
We use [•, •] to denote a column vector, and let 

w(x,t) = [u(x,t),v(x,t)]. 
Let q = (q 1: .., q d ) e SI = (N U {0}) d and let 



e i( x ) = JJcos (g^i) , 



Then {e q (a;)} qg Q forms a basis of the space of functions in T d that 
satisfy Neumann boundary conditions (1.2). We look for a normal 
mode to the linear Keller-Segel system (1.6) and (1.7) of the following 
form: 

(1.8) w (x,t) = r q exp (A q t) e q (x), 

where r q is a vector depending on q. Plugging (1.8) into (1.6)-(1.7) 
yields 

\ r _ ( -W 2 xUq 2 
AqFq ~\ f ~Dq 2 - k 

where q 2 = Ylt^Qi- ^ nontrivial normal mode can be obtained by 
setting 

det(^_T 2 K 7ff +k )^ 
This leads to the following dispersion formula for A q : 

(1.9) A q + {q 2 (fj, + D) + A;}A q + q 2 {fi (Dq 2 + k) - X U f} = 0. 

Thus we deduce the following well-known aggregation (i.e., linear in- 
stability) criterion by requiring there exists a q such that 

(1.10) t,(Dq 2 + k)- X Uf<0, 
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to ensure that (1.9) has at least one positive root A q . This clearly 
implies that /ik — yl) f < 0, and an elementary computation of the 
discriminant yields: 

{q 2 (fi + D) + k} 2 - Aq 2 {fi (Dq 2 + k) - X U f} 

= g 4 (fj, - D) 2 + k 2 + 2q 2 (/i + D) k + 4q 2 {-fik + X Uf} 
> 

for q. Therefore, there exist two distinct real roots for all q to the 
quadratic equation (1.9), which we denote 

A_(q) < A + (q). 

We denote the corresponding (linearly independent) eigenvectors by 
r_(q) and r + (q), such that 

A±(q) + Dq 2 + k 



(1.11) r±(q) 
Clearly, for q large, 



/ 



,1 



H {Dq 2 + k) — X Uf > 0. 

Hence there are only finitely many q such that A + (q) > 0. We therefore 
denote the largest eigenvalue by A max > and define 

^max = {q G such that A+(q) = A max }. 

It is easy to see that there is one q 2 (possibly two) having Aq (q 2 ) = A max 
when we regard A+ as a function of q 2 . We also denote v > to be the 
gap between the A max and the rest. 

Given any initial perturbation w (x, 0), we can expand it as 

w ( x >0) = ^w q e q (a;) = ^{w q r_(q) + w+r + (q)}e q (x), 

q en q Gf2 

so that 

(1.12) w q = w q r_(q) + w+r + (q). 

The unique solution w(x,t) = [u (x,t) ,v (x,t)] to (1.6)-(1.7) is given 
by 

(1.13) w(x,t) = ^{w q r_(q)exp (X~t) +w+r + (q)exp (A+t)}e q (x) 

= e°w (x, 0) . 

For any u(-,t) G [L 2 (T d )] 2 , we denote ||u(-,t)|| = \\u(-,t)\\ L2 . Our 
main result of this section is 
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Lemma 1. Assume the instability criterion (1.10) is valid. Suppose 

w (x, t) = [u (x,t) , v (x,t)] = e ct w (x, 0) 

as in (1.13) is a solution to the linearized KS system (1.6)-(1.7) with 
initial condition w (x, 0). Then there exists a constant C\ > 1 depend- 
ing on k, U, D, fi, f, x, such that 

\\w(-,t)\\ < d exp (\ max t) ||w (-, 0)11 , 

for all t > 0. 

Proof. We first consider the case for t > 1. By analyzing (1.9), for q 
large, we have 

hm = -/i, -D 

q->oc q z 

respectively. Notice that from the quadratic formula for (1.9), 



A+ - A q > 2 y/-fj,k + X Uf 



q 2 



From solving (1.12) 



- ttt — r~\ — ^vrk±(q)l x l w ql 

q det[r_(q),r+(q)] 

< C?|Wq|, 

we deduce that for t > 1 and q large, 

\w^r±(q) exp (Aqt) | < Cg|w q | exp(— min{/i, D}q 2 t) < C|w q |. 

Thus we deduce the Lemma on the linear growth rate for t > 1 by the 
formula (1.13). 

On the other hand, for finite time t < 1, it suffices to derive the stan- 
dard energy estimate in L? . From the Neumann boundary conditions, 
we can take ux (1.6) and add Av x of (1.7) to get 

+ / {/2\Vu\ 2 + AD\Vv\ 2 -xUVvVu} + Ak / \v\ 2 
Ji d Jf d 

Afuv. 

The integrand of the second integral can be chosen non-negative 
(1.14) /i | Vu\ 2 + AD \Vv\ 2 -xUVvVu > ^ \Vu\ 2 + ^ > 0, 



KELLER-SEGEL INSTABILITY 5 

if the constant A is 

(1.15) A = ^l. 

It thus follows that 

and the Gronwall inequality implies 

||w(.,OII <Cexp(C0||w(-,0)||, 
for some C > 0. This immediately implies our lemma when t < 1. □ 

2. Main Result 

Let 9 be a small fixed constant, and A max be the dominant eigenvalue 
which is the maximal growth rate. We also denote the gap between the 
largest growth rate A max and the rest by v > 0. Then for 5 > arbitrary 
small, we define the escape time T s by 

(2.1) 9 = 5exp(\ max T 5 ), 

or equivalently 

'"max " 

Our main theorem is 

Theorem 1. Assume that the set of q 2 = Yl'i=i < l'i satisfying instability 
criterion (1.10) is not empty for given parameters /z, D, k, x, / o,nd U. 
Let 

Wofc) = ^{^ q r_(q) +w+r + (q)}e q (a;). 

G H 2 such that ||w || = 1. Then there exist constants So > 0, C > 0, 
and 9 > 0, depending on k, U, D, ji, /, x, such that for all < 5 < 5o, 
if the initial perturbation of the steady state [U,V] in (1.3) is 

w 5 (x, 0) = 5w , 

then its nonlinear evolution w s (t,x) satisfies 

\\w\t,x) - 5e Amax * ^ w q r +( c e ^( x )\\ 

< C{e~ vt + 5\\w \\ 2 H2 + 5 e W} 5e A max t 

for < t < T 5 , and v > is the gap between A max and the rest of A q 
in (1.9). 
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We notice that for < t < T 5 , <5 e Amax * < 9, is sufficiently small. 
As long as iu+ o ^ for at least one qoGfi max , which is generic for 
perturbations, the corresponding fastest growing modes 

||6e w ™q r +(q)e q || > 5e w |w;q () l|r+(qo)|, 

have the dominant leading order of <5e Amax< . Our theorem implies that 
the dynamics of a general perturbation is characterized by such linear 
dynamics over a long time period of eT s < t < T s , for any e > 0. In 
particular, choose a fixed q G fi max and let 

w {x) = - — - — rre qo (x) 



iqoj 



then if t = T s 



i r +iqoji 



< c{6 v/Xnax +e' 2 }, 



hence 

||w 5 (t, -)|| > - C{5 V,X — + 9 2 } > 9/2 > 0, 

which implies nonlinear instability as 5 — > 0. The instability occurs 
before the possible blow-up time. 

In the early work of Keller and Segel [15] in 1970, they formulated 
the advection-diffusion system (1.1) which consists of two parabolic 
equations and viewed the initiation of Slime mold aggregation as in- 
stability. Linearized system was used to analyze early stage of pat- 
tern formation and its instability around homogeneous steady states. 
This Keller-Segel model has since received much attention and there 
have been many contributions on this subject such as aggregations, 
dynamics of blow-ups, travelling waves. See [1], [2], [3], [9], [10], [11], [8], 
[13], [16], [17], [19], [20], [21] for related results. Linear stability and insta- 
bility of stationary solutions with more general nonlinearity was studied 
in [22] using bifurcation analysis. However, nonlinear evolution of the 
pattern formation has yet been fully understood for the Keller-Segel 
model, to the authors' knowledge. 

We rigorously prove that linear fastest growing modes determine un- 
stable patterns for the full Keller-Segel system (1.4) and (1.5), over a 
time period of the order ln|. Each initial perturbation certainly can 
behaves drastically differently from another, which gives rise to the rich- 
ness of patterns. On the other hand, the dominating linear dynamics 
over a fixed finite dimensional space of maximal growing modes en- 
sures that there is a common characteristic pattern for a general class 
initial data. Therefore, we believe that our result indeed provide a 
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mathematical description for the pattern formation in the Keller-Segel 
model. 

Our paper stems from a program to study various nonlinear instabil- 
ities for non-dissipative systems arising in mathematical physics [5], [6], 
[7], [12], where severe higher order perturbations (unbounded in the L 2 
norms, for instance) occur. Indeed, for many such systems without dis- 
sipation, the passage from linear instability to nonlinear instability is 
very delicate. If there is a dominant eigenvalue, then a bootstrap argu- 
ment was developed by Strauss and the first author to prove nonlinear 
instability, for the perturbation initially along the dominant eigenfunc- 
tion. The key is to try to control the nonlinear growth of higher-order 
energy norm for the perturbation by the linear growth rate, up to 
the time T s . Very recently in [4], based upon a precise linear analysis, 
dynamics of general perturbation can be characterized by the linear dy- 
namics of fastest growing modes for unstable Kirchhoff ellipses. This 
marks a beginning of a quantitative description of instability. 

Our research is inspired by the work [4]. In the presence of dissipa- 
tion, continuum spectra are absent in bounded domain, which leads to 
finite number of dominant growing modes. Moreover, natural higher- 
order energy estimate now can be easily combined with the bootstrap 
idea to control the nonlinear term — %V(itVt> ) in the L 2 space. Since 
our method is general, we believe that such kind of pattern formation 
should exist for a wide class of systems with dissipation. 



3. Bootstrap Lemma 
We state existence of local-in-time solutions for (1.4)-(1.5). 

Lemma 2. (Local existence) For s > 1 (d — 1) and s > 2 (d — 2,3), 
there exist a T > and a constant C depending on k, U, V, D, /i, f, \ 
such that 

\Mt)\\ HS <c \\w(o)\\ Hs . 

We now derive the following energy estimates for <i-dimensional chemo- 
taxis model with d — 1,2,3. 
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Lemma 3. Suppose that [u (x,t) , v (x, t)] is a solution to the full system 
(1.4)-(1.5). Then 

2 dt ^ J T d | 1 D\i 1 1 ' 



\d\=2 



=2 

< Col |w| 1^/2 1 1 V 3 w| |^ + C 2 \\u\ 




where C is the universal constant while C 2 



t/ 6 x 6 / 6 

2D 3 /i 5 k 3 ■ 



Proof. We first notice that the Keller-Segel equation preserves the even- 
ness of the solution w(x,t), i.e., if w(x,t) is a solution, then w(—Xi,t) 
is also a solution. We can regard the Neumann problem as a special 
case with evenness of the periodic problem by standard way of even 
extension w(x, t) with respect to one of the Xi. For this reason we may 
assume periodicity at the boundary of the extended 2T 3 = (— ir, n) d . 
Since now there is no contributions from the boundaries, we can take 
second order <9-derivative of (1.4) and add A x d of (1.5) to get 



[ {/i|V<9w| 2 + AD\Vdv\ 2 - X UVdvVdu} + Ak [ 

J2J d J21 

X [ d{uVv}Vdu +Af f dudv 

J2T d J2T d 



= h + h, 

where the constant A is given in (1.15). As in (1.14), the second inte- 
grand is bounded below by 

^| V ^| 2 + -^|V«%| 2 . 
2 1 1 2/i 1 1 

The nonlinear term l\ is bounded by 

h = J \d{uVv)-Vdu\d* 

< \\u\\ Loo \\Vdv\\ \\Vdu\\ + \\Vu\\ Loo \\dv\\ \\Vdu\\ 
+ \\du\\ l|Vu|| L oo ||V9u||. 
We apply the following the Sobolev imbedding to control ||m|| loo 

(3.1) IMIl°o( 2 T<*) - Co \\g \\ H 2(2T d ) ' 
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for d < 3. Moreover, from the periodic boundary conditions, 

/ Vu = / Vv = 0, 

J2T d J2J d 

we also use the Poincare inequality 

(3.2) \\g\\<\\g\\ L 4 (2T d ) < Co ||V^|| if d < 3, 

to further get 

||Vu|| L oo + ||Vu|| L <» < C {\\Wu\\ H2 + llVull^,} < C ll 9Vw ll> 

\d\=2 

where Co is a universal constant. Hence I\ < Co 1 1 w 1 1 i? 2 1 1 V 3 w 1 1 2 as 
desired. 

Finally, I 2 is simply bounded by 

/A f 2 Ah 
dudv < ^-\\du\\ 2 + — \\dv\\ 2 

By the interpolation between ||V<9w|| and \ \u\ \ , the first term above is 
bounded by 

for any a > 0. We can choose a such that 4*T a = I^ - Collecting terms, 
we conclude the proof. □ 

We are now ready to establish the bootstrap lemma, which controls 
the H 2 growth of w(x,t) in term of its L 2 growth. 

Lemma 4. Suppose that w(x,t) is a solution to the full system (1-4)- 
(1.5) such that for® <t<T, 

ll ( + \\\ <r 1 • fix)' 
\\w(-,t)\\ H 2 < —mm <^ -, 



C 1 4 ' 2/i 

and 

(3.3) ||w(-,0ll <2C ie w ||w(-,0)||, 

then we have for < t < T, 

l|w(t)||^<C 3 {||w(0)||^ + e 2 ^||w(,0)|| 2 } 

where C 3 = C 2 max{%£, ^} x max{g-, 1} > 1. 
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Proof. It suffices to only consider the second-order derivatives of w(x, t). 
From the previous lemma and our assumption for ||w||#2, we deduce 
that for < t < T 

,2 



So that by (3.3) and an integration from to t, we have 
<J2 I l\du(0)\ 2 + ^^\dv(0)\ 2 \ + ^^e 2X ^ *| 



w 



for < t < T. Now our lemma follows directly by separating the cases 
of M > 1 and M < 1. □ 

4. Nonlinear instability and pattern formation 

We now prove our main Theorem 1: 

Proof. Let w s (x, t) be the family of solutions to the Keller-Segel system 
(1.4)-(1.5) with initial data w 5 (x, 0) = 5w . Define T* by 

T* = sup ji | \\w\t) - 5e ct w \\ < ^5exp (A max t) j . 

Note that T* is well defined. We also define 

i J/< (ux)' 



T** = sup It | \\ w \ t )\\ H2 <—mm{^ 



C 14' 2^ 
We recall T 5 in (2.1) where 9 is chosen such that 

Amax At {Ux}' 



(4.1) C C 3 9 < min 



4 ' 8 ' 4/i 



We now derive estimates for H 2 norm of w s (x, t) for < t < 
mm{T*,T s ,T**}. First of all, by the definition of T*, for t < T* and 
Lemma 1 

or* 

\\w s (t)\\ < _i5exp(A max t). 

Moreover, using Lemma 4 and applying a bootstrap argument yields 
(4.2) K(t)|| H2 < v^lKII^ + 5e w }. 
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We now establish a sharper L 2 estimate for w 5 (x, t), for < t < 
min{T**, T S ,T*}. We first apply Duhamel's principle to obtain 

w" 5 (t) = 5e ct w - f e c{t ' T) [V • (u s (r) Vv s (r)) , 0}dr, 
Jo 

Using Lemma 1, (3.1), (3.2), and Lemma 4 yields, for < t < min{T s ,T** 
\\w s (t) -5e ct w \\ 

<d f e Amax( *~ T) || V • (u s (r) Vv s (r)) || dr 
Jo 

<d [ e A — ^ \\u 5 {r)\\ Lao \\V 2 v 5 {r)\\dr 
Jo 

+ d f e w (*~ r ) ||Vu'(r)|| L4 \\Vv s (r)\\ L4 dT 
Jo 

^dCo / e w{ *- r) ||w 5 (r)||^ 2 rfr. 
Jo 

By our choice of t < min{T*, T**, T s }, it is further bounded by 
(4.3) \\w 5 (t) - 5e ct w \\ 

< dCoCg / e x ^ {t -^{5 2 \\w \\ 2 H2 + 5 2 e 2X ^ T }dr 
Jo 

^max '"max 

We now prove by contradiction that for S sufficiently small, 

T s = min{T 5 ,T*,T**}, 

and therefore our theorem follows by further separating q e fi m ax and 
move q ^ il ma x in (1-13) to the right hand side . 

If T** is the smallest, we can let t = T** < T 5 in (4.2) 



||w d (T**)|| H2 < VCs{S\\w \\ m + Se A — 1 } 
= ^{51^011^ + 0} 

for v^Cs^l l w o| |iJ 2 < i^ 111 ^ 11 !!' 21 }' ^ our cn °i ce °f $ i n (4-1) with 
Cq > 1. This is a contradiction to the definition of T**. 
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On the other hand, if T* is the smallest, we let t = T* in (4.3) to 

get 

||w 5 (T*)-5e CT *w || 

< C 1 CoC 3 { ll ^° 11 ^ + J—5e x ™ TS }6e x ™ T * 

-""max -\riax 
^max ^max 

< — be max , 

2 

for C C 3 V H < V 4 f o r small, by our choice of 6 in (4.1). This 

m ax ' 

again contradicts the definition of T* and our theorem follows. □ 
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